Journal of Applied Crystallography

X-ray fiber diffraction patterns exhibit four-fold symmetry that can be exploited, through folding and averaging, to improve signal-to-noise ratio. Accurate folding requires a precise sub-pixel estimate of the symmetry center and precise orientation of the meridional pattern axis to the fiber axis: small center or angular errors blur diffraction features, reduce layer-line sharpness, and introduce errors in spacing measurements. A pixel-level estimate is often too imprecise for this purpose, and detector gaps further complicate the alignment objective. We formulate the masked quadrant-folding problem, define a four-quadrant symmetry loss that consistently excludes invalid pixels, and evaluate several refinement strategies: hierarchical coarse-to-fine grid search; ECC-based rigid registration with global center/orientation correction fitting; ECC registration followed by local gradient refinement; and a hybrid that appends a local grid search on a cropped pattern. Direct gradient optimization from the rough QF alignment was found to be unreliable. Grid search provides a robust, interpretable baseline that directly minimizes the folding objective but is substantially slower than registration; ECC gives a fast near-correct alignment, and the hybrid closes the accuracy gap to brute-force search at a fraction of its runtime. On real datasets with calibration data, applying a calibration center with optimized rotation is effectively optimal. The hybrid center-refinement method has been integrated into the MuscleX package.

Three-dimensional electron diffraction (3DED) is an emerging and useful technique for solving molecular structures of small and biological macro-molecules from nanometre-sized crystals. We present our automated data processing workflow for 3DED datasets collected at Diamond Light Sources electron Bio-Imaging Centre (eBIC). For this purpose, we developed a package called AutoED. The processing pipeline includes data collection, analysis of the beam position, metadata gathering, file conversion, and finally data processing using xia2 (which supports both DIALS and XDS). The processing results are captured in a summary report produced by AutoED. Our main goal is to reduce the workload of electron diffraction scientists, but also to enforce good standards already used in macromolecular crystallography (MX). All the collected 3DED datasets are automatically converted into NeXus data format which is considered a Gold Standard for MX. This standardized data format allows for all the relevant metadata about the experiment to be kept together with diffraction images. We also discuss the methods used in AutoED to determine the electron beam position on diffraction images.

Accurate estimation of the contrast transfer function (CTF) of tilt images is a critical first step in cryo electron tomography (cryoET), enabling reliable recovery of high-resolution structural information from thick, heterogeneous specimens. This challenge is especially acute in in situ cryoET, where macromolecules are imaged in their native cellular environment, often at high tilt and through substantial specimen thickness, with correspondingly low signal-to-noise ratios. Although CTF parameters can be later refined using reference-based approaches, accurate initial estimates are critical for downstream processing and the interpretability of tomographic reconstructions, yet they remain difficult to automate. Here, we present Quinoa, a software package designed to address these challenges. Quinoa first validates the tilt geometry and assesses data quality to generate robust initial estimates of defocus and phase shift. These estimates are then refined through optimization of a single global model, enabling precise fitting of the per-image defoci, tilt-dependent astigmatisms, time-dependent phase shifts, the specimen orientation (rotation, tilt and pitch) and the specimen thickness. Notably, and as a key distinguishing feature of this approach is that Quinoa fits equiphase-binned polar power spectra. This substantially reduces the computational cost of optimization without sacrificing accuracy, enabling more progressive and exhaustive refinement passes that further improve robustness. We validated Quinoa using both simulated and experimental data and benchmarked its performance against Warp, Ctfplotter, CTFMeasure, and AreTomo. Our results show that Quinoa is the most robust approach across all simulated cases, maintaining high accuracy even in the simultaneous presence of severe astigmatism, high specimen inclination and variable phase shift. Integrated recovery mechanisms further allow Quinoa to adapt automatically to a wide range of pixel sizes, defoci, astigmatisms and specimen thicknesses. Despite fitting a more complex and dynamic model, Quinoa remains extremely efficient due to extensive GPU acceleration, making it well suited for real-time monitoring during data collection as well as high-throughput offline batch processing. By improving automated CTF estimation in challenging tomographic data, Quinoa supports more accurate structural analysis of cells and tissues in situ.

Four-dimensional scanning transmission electron microscopy (4D-STEM) enables the acquisition of diffraction patterns at every probe position in a dense array. For imaging applications this approach offers significant benefits in terms of spatial resolution and contrast enhancement. In this work, we present the development of a synchronous scan generator integrated with SerialEM software to enable automation of complex experimental protocols such as tomography. The proposed hardware functions as an interface between SerialEM, the scan controls of the microscope, a fast annular dark-field detector, and a synchronized trigger for a pixelated detector. Our previous implementation, named SavvyScan, relied on a dedicated computer equipped with a multichannel acquisition and signal-generation cards, as well as a separate microcontroller for synchronization. Here, we report a low-cost implementation based on a Red Pitaya board, utilizing direct programming of its embedded FPGA and Linux server components. We provide detailed instructions for system installation and operation, along with practical guidance for modifying the source code. System performance is validated through oscilloscope measurements and imaging of a replica grating sample. The utility of the approach is further demonstrated by generating a 3D electron tomogram of a cryogenic sample of mitochondria from a tilt series of shadow montage projections.

Small-angle X-ray or neutron scattering (SAXS/SANS) analysis of large unilamellar vesicles (LUVs) is often limited by high-dimensional bilayer models and the lack of dedicated, statistically rigorous workflows. Here, we introduce SAS_MoCa, an open-source Python package that integrates a compositional scattering density profile (SDP) description of lipid bilayers with a separated form factor (SFF) treatment of vesicle size and polydispersity, and couples these highly parameterized models to an adaptive thermodynamic simulated annealing algorithm formulated within a constrained Bayesian framework. SAS_MoCa enables users to incorporate quantitative prior information from, e.g., previous SAXS/SANS studies, dynamic light scattering, NMR, or molecular simulations, and returns full posterior parameter distributions, uncertainties (reported as medians and median absolute deviations) and correlations even from single SAXS curves. Validation on POPC, POPE and DMPC SAXS-only data demonstrates that the method yields reproducible structural parameters with uncertainties comparable to joint SAXS/contrast-variation SANS analyses. The modular architecture of SAS_MoCa facilitates extension to additional lipid systems and future joint SAXS/SANS or SANS-only applications.

Chemical imaging is a powerful branch of modern microscopy encumbered by a lack of flexible, high-throughput analysis tools. Bespoke analytical pipelines typically perform ratiometric analysis on two layers in a multispectral image to describe the relative composition of molecules in a sample. This strategy has been implemented across fields, spanning histopathology, cell biology, environmental science, and materials science. The commercialization of chemical imaging microscopes has facilitated the collection of large multispectral datasets, necessitating accessible ways to process them. This paper describes Multispectral Analysis Graphical User Interface (msaGUI), a desktop graphical user interface to analyze individual and batch datasets of multispectral images. Data is loaded as CSV, TSV, or TIFFs and processed through a user-defined sequence of modular image operations that can be flexibly combined, e.g. to reduce spectral crosstalk or background noise. After analysis, data is visualized as exportable images, histograms, and statistics. To yield publication-quality figures, outputted images are fully customizable. Written in Python with open-source libraries, the msaGUI program is packaged into an executable for Windows and Mac for a fully no-code application. Other operating systems are supported via the Python source code. In summary, msaGUI provides a rapid and user-friendly solution for analyzing and visualizing multispectral data.

Directional anisotropy in electron density provides key information about chemical bonding that is not readily accessible from conventional electron-density maps. Here, a model-independent framework is presented for decomposing experimental structure factors into angular components using spherical harmonics. Reciprocal-space projection onto spherical harmonics followed by standard Fourier synthesis yields angularly filtered density maps. The{ell} = 0 component captures the isotropic part of the density, while the{ell} = 1 components resemble px, py and pz-like dipolar functions that highlight directional electronic structure. Applications to high-resolution datasets, including urea, the Gly-Ala dipeptide and a 0.97 [A]{beta}-lactamase structure, reveal chemically interpretable dipolar features associated with carbonyl and amide bonds, N-H interactions and aromatic{pi} systems. Quantitative analysis using bond-centred sampling demonstrates stable dipolar signatures that remain detectable under moderate resolution truncation. These results establish spherical-harmonic angular decomposition as a practical framework for extracting directional electronic information from crystallographic electron-density maps. SynopsisAngular decomposition of experimental structure factors reveals dipolar anisotropy and directional electron-density features that are directly meaningful for chemical interpretation.

Cryo-electron microscopy (cryo-EM) of biological specimens is limited by radiation damage and a low signal-to-noise ratio (SNR). Here, we show that reducing the illuminated area substantially slows the observed diffraction decay in protein microcrystals. We further show that narrow parallel-beam electron diffraction from thin non-crystalline biological specimens provides substantially higher reciprocal-space SNR than conventional cryo-EM imaging. We developed a multimodal scanning workflow, 4D-para-STEM, that records narrow-beam diffraction patterns together with corresponding images. Using viruses, peptide assemblies, and microtubules, we demonstrate interpretable diffraction signals from both crystalline and non-crystalline biological specimens. Together, these results show that narrow parallel-beam scanning reduces observed radiation damage and improves the SNR in cryo-EM.

Cryo-electron tomography is an expanding technology for the study of macromolecules, viruses, and cells. It is often applied to specimens that are too large or heterogeneous for methods based on 2D image averaging such as single particle analysis, e.g., intracellular membranes or organelles. Current practice records a tilt series of projection images in rotation. Reconstruction is normally an ill-posed mathematical problem. Particularly for the under-determined case of sparse data, discrete tilt angles, and a limited tilt range, characteristic artifacts appear in the reconstructed slices. Much of what appears as noise is in fact structural: the projection of contrast from different planes. Various schemes are employed to regularize the reconstruction, including machine-learning frameworks built on neural networks. To the extent that the noise is structural, it might be suppressed by deconvolution with a suitable kernel. This was demonstrated and has been used regularly in cryo-STEM tomography of thick specimens where the under-sampling problem is particularly acute. Here we present 3dcon as an open-source extension of the entropy-regularized deconvolution algorithm that had been adopted from fluorescence microscopy. It takes advantage of modern computing hardware for convenient and fast processing. Deconvolution is entirely algorithmic, meaning that successful processing of the data does not depend on the data itself. As such it should be robust in a wide variety of applications.

Images in the electron microscope are formed by electron scattering and focusing. The spherical geometry of these processes gives rise to two coherent, conjugate spherical wave fronts, known as Ewald spheres. These spheres are associated with the two halves of the contrast transfer function (CTF), and their widths are determined by the focal gradient through the specimen. To properly correct for the CTF, each half of the CTF must be applied to an image individually and integrated into the reconstruction into the corresponding Ewald sphere. Theory indicates that this dual Ewald sphere reconstruction method should recover the maximal amount of information possible. This method was compared to the other reconstruction methods commonly used: the projection approximation (ignoring the Ewald sphere), the simple insertion and the single sideband methods. In simulated reconstructions the dual Ewald sphere method recovered the most information when the correct half of the CTF is matched to the corresponding Ewald sphere. If the wrong half is matched, the result worse than the projection approximation method. Examining reconstructions from real data indicated that the dual Ewald sphere method performs at least as well as the simple insertion method, but not as good as in simulations. The likely reason is the two-fold ambiguity in the assigned orientations of the particle images, which remains an issue to pursue in further studies. In conclusion, the dual Ewald sphere reconstruction method may offer the best way to calculate very high resolution reconstructions when the micrograph quality warrants it. HighlightsO_LIThe dual Ewald sphere reconstruction corrects for the two halves of the CTF. C_LIO_LIThe signs of the two halves of the CTF must correspond to the focal gradient. C_LIO_LIDetermining the focal gradient for individual particle images remains unresolved. C_LIO_LIComplex reconstructions indicate any real space phases are artifacts. C_LI

Accurate assessment of ligand coordinate-density consistency across different resolutions remains challenging in macromolecular crystallography. We introduce the atomic Box Correlation Coefficient (aBCC), an atom-level metric for evaluating the consistency between ligand atomic coordinates and electron density in a resolution-standardized framework. To predict aBCC values from electron-density maps, we developed QAEmap, a machine-learning model based on three-dimensional convolutional neural networks (3D-CNNs). The model was trained using Fourier-truncated electron-density maps and corresponding ligand coordinates generated from high-resolution structures in the Protein Data Bank. It was evaluated using both Fourier-truncated electron-density maps and experimentally determined PDB structures. was evaluated using both Fourier-truncated electron-density maps and experimentally determined PDB structures.The prediction accuracy gradually decreased with decreasing resolution, but remained reliable up to [~]3.5 [A]. These results demonstrate that aBCC enables resolution-standardized atom-wise evaluation of coordinate-density consistency across different resolutions and provide a foundation for further development and refinement of machine learning-based coordinate validation. SynopsisWe introduce the atomic box correlation coefficient (aBCC), a machine learning-based metric for the resolution-standardized atom-level evaluation of ligand coordinate-density consistency in crystallographic structures. aBCC provides a common framework for assessing and communicating the local coordinate reliability between structural biologists and researchers in structure-based drug discovery.

Damage to biological specimens by the electron beam is the fundamental resolution-limiting factor in cryoelectron microscopy (cryo-EM) single particle analysis. There is, however, currently no method to accurately infer fluence-dependent changes to the specimen structure during electron irradiation. We develop a Bayesian framework to fit a sequence of atomic models to a series of cryo-EM reconstructions produced at increasing fluence. In particular, our algorithm is able to infer the ensemble average position and atomic displacement parameter of every atom in the macromolecule as a function of fluence. Application of the algorithm to cryo-EM datasets shows that the molecule expands during imaging and identifies environment-dependent variations in beam-induced damage. We use our results to propose a stochastic process model of this phenomenon. We envisage that our method will lead to a better mechanistic understanding of radiation damage to biological specimens and may contribute to efforts to mitigate its effects.

This paper presents and critically reviews the results of original and some literature based experimental studies conducted by the authors last years on the structural organization of DNA in dormant (starvation stress), anabiotic dormant (4 HR treatment) E. coli cells, as well as the K12 {Delta}dps strain, which lacks the Dps protein (Dps null E. coli). The experimental data includes small-angle synchrotron radiation diffraction (SAXS) and transmission electron microscopy (TEM) data. Synchrotron radiation diffraction experiments on K12{Delta}dps cells allowed us to conclude that peaks at 44.3, 22.1, and 14.8 angstrom resolutions are associated exclusively with ordered DNA organization. Peaks at 44.3, 22.1, and 14.8 angstrom resolutions are also observed for samples of dormant (starvation stress) cells and anabiotically dormant cells. Therefore, this ordered DNA organization also applies to samples of dormant and anabiotically dormant cells. A model is proposed that considers the ordered DNA organization in the cell as a cholesteric liquid crystal. The powder diffraction pattern calculated based on this model is compared with experimental small angle X ray scattering (SAXS) data obtained on Dps-null cell samples. The model completely reproduces the key features of the experimental diffraction pattern from Dps-null cell samples. Accordingly, the cholesteric liquid crystal model corresponds to DNA packaging in dormant and anabiotically dormant cells. Cholesteric liquid crystal ordering should be further considered in all models of cellular DNA packaging. To address the question of which structural organization of DNA predominates in the cell: the cholesteric liquid crystal or nanocrystalline or whether they coexist and fully manifest themselves under different external conditions, it is necessary to utilize the latest methodological advances in structural analysis.

Time-resolved serial crystallography (TR-SX) has emerged as a powerful method for capturing ultrafast structural dynamics in proteins. TR-SX continues to produce remarkable studies, revealing previously unobserved transient states and providing deeper insights into processes such as drug targeting, DNA repair, and photosynthesis. However, extracting weak structural signals from noisy time-resolved datasets remains a major challenge. Robust computational methods are therefore required to isolate the signals associated with the underlying transient states. Importantly, this should be performed in reciprocal space to preserve compatibility with established downstream structure refinement workflows. Here, we introduce a framework for kinetic decomposition directly in reciprocal space that enables separation of kinetically distinct structural states. The method decomposes crystallographic data according to a predefined kinetic model, improving the recovery of weak transient signals and enhancing mechanistic interpretation from limited time-resolved datasets. We validate the framework using simulated data based on a previously published time-resolved crystallography study and demonstrate its application to a new TR-SX dataset comprising 17 time points. We show that the method separates the reciprocal space signatures of four intermediates by incorporating kinetic information from a predefined reaction model. This establishes a workflow for extracting kinetic states directly from time-resolved X-ray diffraction data that can be seamlessly integrated into existing crystallographic structure-determination pipelines.

Microscopes are essential tools for discoveries on a scale invisible to the unaided human eye. The development of immuno-fluorescence followed by molecular biology techniques and fluorescent fusion proteins have revolutionised the use of optical microscopy in bioscience. The quality of the data produced is dependent upon the sample, its preparation and the instrument used. However, instruments can degrade over time without easily visible changes to the produced images and, in turn, negatively impacts results. By testing instruments and doing comparisons between results over time and between different instruments, problems can be highlighted and corrective action can be taken. Using small fluorescent beads the point spread function (PSF) of the microscope can be recorded and the image resolution measured. Beads were prepared in a concentration matched to the field of view size and dried onto coverslips and mounted on slides. The beads were then imaged as 3D Z-stacks of sufficient size to fully enclose the PSF of the system. This data was uploaded to OMERO and processed using OMERO-metrics, an OMERO plugin developed for this purpose. This paper summarizes the development of workflows and protocols to enable this process, presents the results obtained and demonstrates the detection of significant instrument issues.

Drug discovery remains one of the costliest and most time-intensive endeavors in the pharmaceutical pipeline, with average development costs exceeding $2.3 billion per drug, timelines spanning more than a decade, and attrition rates above 90% in clinical trials. While computational methods have expanded the searchable chemical space, current pipelines remain fragmented and largely inaccessible to researchers without deep interdisciplinary expertise. Here we present F.A.D.E. (Fully Agentic Drug Engine), a multi-agent, open-source platform that converts natural language queries into potential drug candidates, substantially lowering the expertise barrier to advanced computational drug discovery. F.A.D.E. employs a three-branch hierarchical architecture that adapts to the level of available structural data for any protein target, integrating structure prediction, binding pocket detection, equivariant diffusion-based de novo ligand generation, and binding affinity estimation into a single automated pipeline. We validate F.A.D.E. on two structurally distinct targets: the epidermal growth factor receptor kinase domain (EGFR), a well-established oncology target, and cellular retinol-binding protein 1 (CRBP1), a lipid-binding protein involved in retinoid metabolism. For EGFR, our generated candidates achieved QED scores of 0.85 compared to 0.46 for the co-crystallised reference ligand, demonstrating marked improvement in predicted drug-likeness. Results across both targets confirm that F.A.D.E. can reliably generate chemically tractable, drug-like hit compounds across diverse protein classes from simple natural language input.

Here, we present MC-Bayes, a Python-based script for processing cryo-electron microscopy EER movies on one or more GPUs using MotionCor3 in a user-friendly manner. Further, it generates the .star files necessary for RELION to perform Bayesian polishing (a.k.a.: reference-based motion correction) with EER movies. Until now, Bayesian polishing of EER data was only possible if the CPU-based "RELIONCor" implementation of MotionCor2 was used, which is sub-optimal on GPU-heavy cryo-EM processing systems. This wrapper was created for those facilities and/or users who may have (many) powerful GPUs, but for whatever reason have few CPU cores or less system RAM. Leveraging MotionCor3, MC-Bayes allows motion correction of EER data 2 or more times faster (depending on system) than the RELION CPU implementation, except in circumstances where dozens or hundreds of CPU cores with high quantities of system RAM can be utilised.

Reported relative biological effectiveness (RBE) values for low-energy X-rays disagree, assays scoring initial DNA double-strand breaks (DSBs) returning about 1.1 and chromosome-level assays 2 to 4. Whether radiation quality varies within the diagnostic range, and how its comparison with a megavoltage reference depends on target scale, has not been quantified on a tube-potential series. A tungsten-anode tube with 1 mm Be and 2.5 mm Al filtration, with copper added in some cases, was modelled in PHITS for 40 to 200 kV. The spectra were transported into a water phantom in which absorbed dose, lineal-energy densities and cluster size distributions were scored for target diameters of 3 nm to 1 micrometre against a cobalt-60 reference; DSB yields were computed in the electron track-structure mode with the PHITS DNA damage tally. Between 40 and 120 kV the depth-dose ratio changed by a factor of 5.7 and the tube output by a factor of 42, whereas the dose-mean lineal energy varied by 2.5 % at 1 micrometre and 1.2 % at 3 nm against a reproducibility of 0.3 %. Relative to cobalt-60 it was 2.05 times larger at 1 micrometre but only 1.08 times larger at 3 nm, while DSB yields per unit dose were 5 to 7 % higher and constant across the range within the 2 % bound set by the statistics. Tube potential therefore changes the amount and distribution of dose but not its physical quality, and a stated RBE is incomplete without the target scale implied by the endpoint.

We demonstrate the use of fluorescence correlation spectroscopy (FCS) to characterize fluorescently-labeled protein production. We use cell-free protein synthesis to express the protein YFP-CopB, a fusion of Chlamydia Outer Protein (Cop) B and Yellow Fluorescent Protein (YFP). CopB is a [~]50 kDa protein believed to have a critical role in chlamydial infection.1 After adding a plasmid encoding YFP-CopB to an E. coli cell-free lysate, protein expression begins. We track the cell-free reaction over several hours using the EI-FLEX, a commercial instrument with FCS capability. As protein is expressed over time, YFP-CopB increases in concentration, and the EI-FLEX detects an increase in fluorescent signal above the background of the cell-free lysate. The FCS data collected gives information about the size, aggregation tendencies, rates of production and fluorescent protein maturation, and concentration of the YFP-CopB produced. The use of FCS concurrent with cell-free synthesis presents a simple method to characterize proteins of interest as they are produced without the need for purification.

Cryogenic electron microscopy (cryoEM) is now routinely used for high-resolution structure determination of biological macromolecules. However, many biological specimens exhibit varying degrees of preferred orientation on cryoEM grids, resulting in uneven sampling of three-dimensional Fourier space. This orientation bias produces anisotropic reconstruction artifacts and, in severe cases, can exacerbate particle misalignment during iterative refinement, thereby limiting the success rate of near-atomic resolution cryoEM structure determination. This protocol provides a practical guide for applying spIsoNet, a self-supervised deep-learning method, to mitigate preferred-orientation issues in cryoEM reconstructions. We describe two complementary workflows: (1) map Anisotropy Correction to correct anisotropic artifacts of cryoEM maps and (2) particle Misalignment Correction, which integrates spIsoNet with RELION external reconstruction to improve particle-pose estimation. We demonstrate these workflows using two influenza hemagglutinin (HA) trimer datasets representing moderate and severe degrees of preferred-orientation bias. The protocol includes installation instructions, parameter-selection guidance, quality-control checkpoints and troubleshooting advice, and can typically be completed in ~7 hours on a workstation equipped with four NVIDIA A100 GPUs. Together, these workflows provide step-by-step guidance for using the open-source spIsoNet software to mitigate the preferred-orientation problem directly from experimental data.